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Abstract 



Bound-state corrections to J/tp production from almost real photons 
are calculated in the colour-singlet model. A systematic, gauge- invariant, 
theory of hard quarkonium processes is used upto 0(v 2 ), where v is the 
relative velocity of the quarks. The internal structure of the meson is char- 
acterised by two parameters, e^/M and V 2 (p(0) / M 2 (p(0) , in addition to the 
usual wavefunction at the origin 0(0). These parameters are constrained 
to be consistent with measured leptonic decay of the J/ip and hadronic 
and radiative decays of r\ c . The calculated corrections to the colour-singlet 
model, which include radiative effects, improve agreement with the exper- 
imental data. 



2 



The main production mechanism for J ftp particles in the inelastic process 
7 + g — > J/ip + X is believed to be the fusion, 7 + g — > J /ip + g. Because of this, 
J I ip production is an important tool for exploring the gluonic distribution inside 
nucleons. The first calculation of this in the so-called "colour-singlet" model 
was performed by Berger and Jones over 15 years ago. Subsequently, several 
authors applied the model to data as they became available. Recently radiative 
corrections to the basic model were calculated by Kramer et al.0, who found 
these to be large at moderate photon energies E 1 ~ lOOGeV. Earlier, relativistic 
corrections had been estimated by Jung et al.H using a model proposed by Keung 
and Muzinich 0. These authors found the corrections to be fairly substantial, 
especially in the high-z (z > 0.8) region where the validity of the colour-singlet 
model is suspect. Inelastic photoproduction of the J/ip has been reviewed by 
Alii- 

The purpose of this letter is to explore the effect of the binding of the 
quarks upon J/ip photoproduction in the colour- singlet model. In the original 
calculations PJ this was totally neglected and the c, c were put on their mass shells. 
This is a sensible starting point because the binding energy e# = 2m — M, and 
the quark relative velocity v, are small parameters: e#/M <C 1 and v 2 /c 2 <C 1. 
Since the c quark mass is only ~ 1.5 GeV, substantial corrections could exist. 
This is indeed suggested by the calculations of refs.[§, §]. However, we do not 
find the calculational method convincing for two important reasons. First, the 
model of ref|| does not treat gauge-invariance satisfactorily. Second, it incorpo- 
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rates binding energy corrections but not wavefunction corrections. Additionally, 
systematic improvement of the model seems difficult. Therefore a re-examination 
of bound-state corrections is important. 

We have recently developed a systematically improvable gauge-invariant 
formalism for the one and two photon (gluon) decays of heavy quarkonia |j) . We 
extend and apply this here to 7 + g — ► J /if) + g. This involves three bosons and 
is therefore considerably more complicated in computational terms. As noted in 
M, the J /ip internal structure is described by more parameters than simply (f)(0), 
the quark wavefunction at zero separation. These are e#/M and V 2 0(O)/M 2 0(O). 
Our starting point is that the photoproduction amplitude 7 + p — > J /if) + X is 
given by the sum of all distinct Feynman diagrams leading from the initial to the 
final state (Fig. la). Each diagram can be written as an integral over the loop 
momenta which, for the lowest order diagram illustrated in Fig. lb, is 

T ™ = / Wf Tr W)^™^)] • (!) 
The tensor f/" AtlM2M3 (/c) is the amplitude to produce a free gluon and two quarks, 
not necessarily on their mass-shells, from a photon and gluon. We call this the 
"hard" or perturbative part, and its expression can be read off from Fig. lb and 
permutations. The "soft" part is the zero-gluon, non gauge-invariant, Bethe- 
Salpeter amplitude, 

M{k)= J d 4 xe lk - x (0\T[if)(-x/2)if)(x/2)}\P,e). (2) 

In equations 1-2, x M is the relative distance between quarks and is the relative 

4 



momentum. 

The next category of diagrams contain a single gluon exchange between the 
blob and one of the two hard propagators (Fig. 2). These all have the general 
form, 

7;™ = / 0ji Tr \ MP ^ mr mi ^ (3) 

Again H PlfX2Pa may be directly read off from the diagrams, and the soft part is, 

M p (k) = J d A xd A ze ik - x e ik '- z ^\T[i){-x/2)A p {z)i){x/2)}\P,e). (4) 

This is a matrix in colour space since, A p =^\ a A ap . The soft gluon which originates 
from the blob has its momentum k' bounded by R~ l ~ fc' « M, where R is the 
meson's spatial size. The two gluon diagram can be included in the same way, 

T 2 ™ = J ^^^TrM^'ik, k', k")H^(k, k\ k"), (5) 

with, 

M p , p ,{k,k\k") =Jd 4 xd*x'd*x"e i{k - x+k ' x ' +k ''^ e). 

(6) 

The gluon self-interaction diagram in Fig. 3b is similarly included. 

As the next step, the hard parts are expanded in the quark and gluon 
relative momenta. The various amplitudes are combined and the Ward identity 
d°S F = —S F ^ a S F is freely used. The upshot of the calculations is that the 
ordinary derivatives combine with gauge fields to yield covariant derivatives, i.e., 
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a gauge- invariant result for the amplitude for 7 + g — ► J/ip + X is obtained, 
(T + Ti + T 2 ) w ^ 3 = Tr[(0|^|P,e)^^ 2 ^ 

+ (0$iB a »'B/3 e)-<9 c *<9 / V 1 '^ 3 

+ (0|^fy|P, e>~^##™ + ...]. (7) 

In the above, h^ 1 ^ 2 ^ = iP* 1 ^ 3 (& = 0). The last term shall not concern us here 
since it is of higher order than v 2 . 

To proceed, one can perform a Lorentz and CPT invariant decomposition 
of each of the hadronic matrix elements in Eq.[7[ This is somewhat complicated 
[0] and involves a large number of constants which characterize the hadron. Con- 
siderable simplification results from choosing the Coulomb gauge, together with 
the counting rules of Lepage et al.0. The result of using this analysis is that, in 
this particular gauge, the gluons contribute at 0(v 3 ) to the reaction 3 Si — ► 7 + X 
and hence can be ignored. Even this leaves us with too many parameters, and 
forces us to search for a dynamical theory describing the essential dynamics of a 
QQ system going beyond the usual non-relativistic potential models. A possible, 
but by no means unique, description is provided by the Bethe-Salpeter equation 
with an instantaneous kernel. This has been conveniently reviewed by Keung 
and Muzinich^] and we shall use their expression for the B-S amplitude in terms 
of the non-relativistic wavefunction <p(p). By projecting appropriately from their 
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wavefunction it is readily established that for 1 states, 



3M 2 1 M 



Using Eqs.[7| and |[ the differential crossection for the subprocess 7 + g — > 
J/ijj + g comes out to be 1 , 
da 256 

— = — a e a 2 s 7i 2 Me 2 q \(f)(0)\ 2 [r]ofo(s,t,u) + rj B f B (s,t,u) + i] W f w (s,t,u)]. (9) 

In the above s, t, and u are the partonic level Mandelstam variables (we omit the 
usual carets) which obey the relation s + t + u = M 2 . The appropriate average 
over initial gluon colours and sum over final gluon colours has been made. If 
radiative corrections are ignored, 

Vo = 1,Vb = w Vw = w (10) 

The function f is the standard, leading order, result: 

, _ (s 2 t 2 + t 2 u 2 + u 2 s 2 + M 2 s t u) 

f {s, t, u) - ^ _ _ M2)2 ^ _ M2)2 . (11) 

The binding energy and wavefunction corrections, fs and fw respectively, are 
slightly more complicated: 
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fB(s,t,u) = [D 



-7 st u(s" + t 4 + u 1 ) + 7M 2 (s 3 t 3 + t 3 u 3 + u 3 s 3 " 



1 We used Mathematica ||, supplemented by the HIP package |]To|, for computation of traces 
and simplification of algebra. 



+ (s 2 t 2 + t V + u 2 s 2 ) (s 3 + t 3 + u 3 + 15s t u) 
+ M 2 stu(s 3 + t 3 + u 3 ) + 29M 2 s 2 t 2 u 2 } , 



(12) 



and, 

fw(s,t,u) = -^-\uistu{s 4 + t 4 + u 4 ) -85M 2 {s 3 t 3 + t 3 u 3 + u 3 s 3 ) 

- 27{s 2 t 2 + t 2 u 2 + u 2 s 2 ) {s 3 + t 3 + u 3 + ^-stu) 

- 139 M 2 s t u(s 3 + t 3 + u 3 ) -463M 2 s 2 t 2 u 2 ]. (13) 

The denominator .D is^, 

D = (s - M 2 ) 3 (t - M 2 ) 3 (u - M 2 ) 3 . (14) 

The total 7p crossection is obtained by convoluting da/dt with the gluon distri- 
bution G{x) in the proton. 

^ G{x)% (15) 



dxdt dt 
Integration of Eq.[| over t in the interval M 2 — s to yields, 

cr(s) = ^--K 2 a e a 2 s e 2 q ^-j^r [r] F + r] B F B + r] W F w + rad.cor.}, (16) 



2 We note that Eq.^ reduces to Eq.23 of Jung et al.|| if the condition rjw — \t\b is imposed. 
This latter condition is equivalent to jjW 2 c/)(0) — ^es<P{0), which is the Schrodinger equation 
for quark relative motion in a potential which vanishes at zero separation. It is also worthy 
of note that the same condition emerges as a renormalization condition in the treatment of 
positronium by Labelle et al|llj] (see their Eqs.ll and 12). However, in our treatment there is 
no principle which apriori constrains t\b to bear a fixed relation to rjyy and therefore both will 
be considered adjustable parameters. 
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where, 



F = 
Fs = I 



l-Ai + 2e + e + 2f - 2f (1 + 2£ + hi 1 ) log^j /(l - 0T*(1 + 
-2 + 16f - 10^ 2 + 48£ 3 + 10£ 4 - 64£ 5 + 2£ 6 

i-3e + we 2 - ioee 3 + 17£ 4 - 5ie 5 ) log el /(i - 3 e 2 (i + 4 , 



and 

F W = \ 



26 - 14f + 210£ 2 - 134£ 3 - 274£ 4 + 150£ 5 + 38£ 6 - 2£ 7 



+ (27+50^ + 257r-292^ + 205r-78^ & -4ir)log^/(l-0 r(l+0 • ( 17 ) 

In the above, £ = s/M 2 . The curves for Fq, Fb and i*V as a function of £ are 
shown in Fig [4]. Note that Fw is a large negative number for very small values 
of s which rises and becomes positive for y/s > 4.75 GeV while Fb is a negative 
quantity for all values of the incoming photon energies. The radiative correction 
to er(s) has already been calculated by Kramer et al.||. They have taken into 
account the modification of the initial gluon densities as well as the box diagrams 
and the splitting of the final gluon into gluon and light quark-antiquark pairs. 
Like the behaviour of the scaling functions plotted in 0, the corrections presented 
in Fig. [4] are also large at moderate photon energies but decrease with incoming 
energies. 

The differential crossection da/dz calculated for y/s = 14.7GeV is shown 
in Fig. [5] where z = Ej/^/E^. We use a simple gluon distribution function 
xG(x) = 3(1 — x) 5 . For the numerical evaluations we take a s = 0.19 and m = 
lA3GeV. Additionally, we take rj w = -.073 and \Rj/^\ 2 = 0.978GeV 3 . This set 
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of parameters is the same as given in Ref. 0. These values, when inputted into 
the theoretical formulae, with radiative corrections evaluated at /x = m yield the 
following values for the decay widths @, 

r(J/V> — >e + e-) = 5.61 KeV 
T(r/ C — ► hadrons) = 9.99 MeV 

T(r] c — ► 2j) = 6.48 KeV (18) 



which are quite close to the values of the experimentally measured decay widths [|T5[] 
5.36 ± .28 KeV, 10.3 ± 3.6 KeV and 8.1 ± 2 KeV respectively. As remarked in 
Ref . , the value of a s chosen for the numerical calculations differs from its value 
deduced from deep inelastic scattering, a s {m c ) w 0.3. Larger values inserted into 
the expressions for the decay rates result in large differences between the wave- 
functions at the origin of J /if) and r] c , violating the assumption that these are of 
0(v 2 ). 

In Fig. [5], we show various theoretical results for da/dz, together with the 



experimental data taken from the EMC and NMC|12|, [13], Ti\. The dash-dot line 
is the zeroth order term calculated by setting r\ B = r\ w = in Eq.||. da/dz 
calculated from Eq.H with the identification r)w = ^Vb, but without radiative 
corrections, is shown with a dotted line (see footnote 2). The dashed line is the 
differential crossection for rjw = —.073 without radiative corrections. Finally the 
solid line is the full result (including radiative corrections) calculated upto 0(v 2 ). 
A K factor of 3.5 has been used, obtained by fitting the crossection including 

10 



all the corrections upto 0(v 2 ) to the experimental dataP[ The overall effect of 
including these terms is that they describe the experimental data quite well in 
the region 0.5 < z < 0.9. 

In Fig.[6a]-[6f], the double differential crossection da/dzdp 2 is shown for 



different z-bins and a comparison is made with the experimental data [|12|, [13 
We see that the relativistic corrections move the theoretical predictions in the 
right direction. However the theoretical curve does not describe the data for the 
highest z region (0.95 < z < 1) where the elastic J/ip production is expected 
to be important. Our formalism works best for the moderate z region i.e, 0.7 < 
z < 0.95. Note that the validity of the colour singlet model hinges essentially 
upon all propagators being hard, a condition which only holds away from the end 
points. Gluonic vacuum fluctuation, a non-perturbative effect, may otherwise be 
important |I6|. 

In conclusion, we have estimated the 0(v 2 ) corrections to the inelastic J/ip 
production from 'jp collisions and shown that the agreement of theory with ex- 
periment is improved if these corrections are included. 



3 The value of the K factor depends upon the value of a s . For a s = 0.3, it is lowered to 1.25. 



11 



Acknowledgements 

We thank D.Wyler for identifying an important error in the preliminary ver- 
sion of this paper, and for numerous comments. This work was supported in part 
by funds provided by NSF Grant No. INT-9122027. H.K gratefully acknowledges 
financial support from the Akhter Ali Memorial Fund for her doctoral research. 



12 



Figure Captions 



Figure: 1 a) The amplitude 7 + p — > J/ip + X, b) One of the 6 leading order 
diagrams. 

Figure: 2 a) One gluon exchange diagram. 

Figure: 3 a) Two gluon exchange diagram, b) Three gluon vertex diagram. 
Figure:4 The functions F Q , Fb and Fw as a function of £ = s/M 2 . Note that 
their contribution is large at moderate incoming photon energies and decreases 



as s increases. 



Figure: 5 The differential crossection da/dz plotted as a function of z at -^i = 
14.7GeV. The solid line is the full crossection (including the radiative correction 
calculated by Kramer et al.[[|). da/dz calculated from Eq.[9] with rjs = —0.076 
and r]w = —.073 is shown by the dashed line while the dotted line corresponds to 
the curve with the choice r] W = \r\ B (see footnote 2). The curve represented by 
dash-dot line is the crossection with r) B = and rjw = 0. A K factor of 3.5 has 
been used to account for the overall normalization. The data points are taken 
from EMC and NMCp], H 0. 



Figure:6 The double differential crossection da/dz/ 'dp 2 for different z-bins. The 
solid curve is predicted from the model including t/b and rjw corrections (K = 
3.5). The dotted curve is for the choice r/w = \f]B- The dashed curve is for 
T] B = and rj W = 0. 
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